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THE PROBLEM OF THE PROPELLER IIT YAW 
WITH SPECIAL REFERENCE TO AIRPLANE STABILITY* 
By Franz Misztal 

1. CALCULATION OF THE FORCES PRODUCED 
BY A PROPELLER IN YAW 



Assume the propeller slopes at angle cp to the hori- 
zontal (fig. l) while exposed to a horizontal air stream 
of velocity w. Then the velocity w can he divided into 
two components: one axial, perpendicular to the plane of 
the blade and one transverse, (in the vertical plane) per- 
pendicular to the axis of rotation. Now, the rotation of 
the "blade ahout the axis at angular velocity oo, "brings 
the true relative peripheral velocity of a "blade element 
into the position {r.d-) r oo + u^ = ro) + u^^ sin d- 

(fig, 2), which, geometrically com"bined with the axial com- 
ponent, yields its relative velocity of flow- The flow 
ahout a hlade element, which is variable across the periph- 
ery, is further disturbed by the self-induction and the 
problem is to define the induced additional velocities and 
through it the true flow attitude about the element in ev- 
ery point of the plane of the blade wherein the forces, set 
up by the element, are given according to the airfoil 
theory. 

A first attempt to explore the energy conditions of a 
propeller in yaw was made by Clarke. (Reference 1.) Dis- 
regarding the induced additional velocities he educed a 
force perpendicular to fiie axis of rotation which was pro- 
portional to the angle of setting of the propeller axis 
with respect to the direction of flight. 



*"Zur Frage der schrag angeblasenen Propeller." Abhand- 
lungen aus dem Aerodynamischen Institut an der Technischen 
Hochschule Aachen, No. 11, 1932, pp. 5-27. 
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Hore recently the problem was attacl^.ed "by s. Fistolesi 
(reference 2) in his invest i.,at ion on a propeller of con- 
stant blade width and constant ^ceomet r i cal coefficient of 
advance on the premises of co:istant induced axial and ve- 
rir-heral additional velocities across the entire blade area 
and distribution of the inducer! axial additive velocities 
caused by the yaw proportional to the radius ?nd confor:.!- 
ably to a sine curve on the periphery, while these become 
zero in the peripheral direction. Owing to the somerhat 
arbitrary mathematical assumption of this analysis which 
probably is to be valid for any yaw of the propeller afvainst 
the direction of the relative air, the discrepancies, at 
least at small angles of attack, that is, normal service 
conditions of the propeller, are - as proved elsewhere - 
abnormal with respect to the actual flow attitude. 

On the other hand, by virtue of the ^.eneral i zat i Oii of 
the latest results of the propeller theory, the flov.^ about 
the propeller blader, is now amenable to exact analytical 
definition. 

The basis of the coordinate system is formed by three 
axes perpendicular to each other (figs. 1 and 2); axis y 
of the propeller, (rotc?tion), axis x in the plane of the 
blade (rotation), perpendicular to the direction of air 
flow a.-'d axis z perpendicular to x. In t'lis reference 
system the velocity w is parallel to the yz plane and 
forms with axis y the angle cf » , which is eqtJivalent to 
the slope of the propeller axis against the riirection of 
air flow and which is called Wangle of attach of the propel- 
ler." 

JIow consider the flow attitude about blade element 
t dr at distance r from the axis of rotation (t = blade 
width). The inflow velocity (or relative fli,(?,ht speed) w 
is divided into axial and transverse components Wq and 
u (fig. l) whereby the axial velocity w^ is held con- 

^0 

stant while 9, "the angle of attack of the propeller" is 

changed. The lateral velocity (sideslir), u yields a 

^o 

tangential component Uq (fi^- 2) constant across the ra- 
dius, and a radial component which, being small compared to 

the axial velocity, can be disregarded. Let hi and ~~ 

represent the induced axial and peripheral additive veloc- 
ities at (r, ^3.) of the plane of the blade, whose result- 



IT.A.C.A. Technical Memorandum No. 696 



3 



ant is — . Then the flow attitude is defined hy the ve- 
2 

locity components , Uq r oo and xi^ and the induced 

additive velocity ^— . (Fig. 3.) 

2 

If the axis of rotation is parallel to the direction 
of flow, that is, Uq = 0 and w - Wq , the corresponding 

induced additive velocities would "be — — , — ^ and — 

2 2 2 

S imul t aneo^isly it is assumed in accord with the g'eneral 
propeller theory that the induced velocities at infinity "be- 
hind the plane of the "blade attain to twice the value even 
in the case of propeller in yaw, or in other words, that 
the small frictional losses on the "blades and the jet con- 
traction and autorotation of the propeller slipstream are 
negl ect ed . 

In agreement with the airfoil theory (fig. 3) the 
thrust and the tangential force of a "blade element at the 
point (r, t9-), relative to area rd^^dr, are: 

dS = [dA cos (3^ - dW sin p» J M = 



2tt 

ztp sin - J drd^ 

2n ^ 2 ^ [tan p» J 

dT - [dA sin + dW cos ] ~ = 



(1) 



2n 

2 



-- Ztp ^- Sin p' [l + 
2tt 2 L tan 3 1 



drd^^ (2) 



where zt = total "blade width at radius r, v = effective 
air flow velocity, = induced pitch angle, c^^ - lift 

coefficient and e = lift-drae ratio (for infinite span) of 
the element , 



According to Figure 3; 

sin r^' = 



^o_+_2 

V 



and 



w. 



. ^ , _ _o_+ 2 1 

2 w„ + w' 1 w^ + wi 

° 2" 2 
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Put t ing ; 



and 



yields 



a 



^0 -S"- 



T7 + ^^0 

0 2 



11 = w tan cp sin S- - Wq^, 



s in ^ ' 



_ 1 



(3) 
(4) 

(5) 



where in 



tan p' = 



X 



Ti ' i X + Ti ' Q_ >v}j: 



(3) 



. = 1 U 



outside radius) an(3 \ = = coef- 



ficient of propeller advance referred to the axial inflow? 
velocity and v is approximately replaced by the undis- 
turbed relative velocity c (assuming — small with re- 

spect to c (fig. 3); Tj ' ^ and r| * . correspond to the 
theoretical and the total induced iciency* of the blade 
element in the propeller theory. 



The additive tangential velocity u^ 
fective angle of air flow of the elem.ent 
Ai • Accordingly its lift coefficient is 

dc 

c^ = c^ + Ai --■rr^ 
a a^ di 



changes the ef- 
i by the amount 



Expr ess ing 



di 



- p ( i in radians) 



* The formulas and concepts used, and presumed as knovrn, are 
from Bienen and v. Karman * s report: On the Theory of Pro- 
pellers, V.D.I., 1924, p. 1237. 
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gives 



[1 + P Ai] (7) 

0 



Making due allowance for the last relations in equa- 



tions (l) and (2) we have, with 
c^ 



w ^ 

p = q : 



dS = 



(1 + p Ai) Eq^zt 



w, 



dT = (1 + p Ai) Rq^zt 



2TTri ' 



a 



k 

1 + e 



X + ri' k - e 

' a 



dxd ^ , 
(!•) 



i X + ri' k^|dxd-5 
A. a /'J 

(2') 



which leaves the change in angle of air flovj Ai relative 
to cp , relative to the position of the "blade element in the 
plane of the propeller and relative to the characteristic of 
the hlade profile. According to figure 3: 



(8) 



or, when replacin;:; the small angle differences by sine 
values: w 



^q 



Ai 



2 c, 



-il- + 



(8') 



The induced velocity c' can he defined "by the momen- 
tum theory, wherehy the small frictional losses, the Jet 
contraction and the autorotation of the vortex system are 
disregarded. The velocity in point {v yd-) is visualized as 
symmetrically distributed on the periphery and we v/rite: 

dA =: c^ (1 + p Ai) ztdrp ^~ = p2Trrdr(w + ~)c^ 



and for the case of the propeller in axial flow of velocity 
Wq v;ith equal angular velocity o) : 



w ' 



dA. 



..^ztdrp _i 



pSirrdr (w^ + ""p"^ 



Again v and v^ are replaced "by c and c^, 
that these two equations finally "become 



so 
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= (1 + pAi) V ^' 
^ 0 c 0 



wix ere 



0 



w„ + — 

t = 



° 2 



2 

Then, eqtiations (8) and (9) concede; 



^ i_ip_ 2_coJ^Co 



1 + p £1° JL t 

2 Co Co 



According to Figure 3: 



-0 ^ ■• 0 
ar-d, "by approximation 



2 . S 



c- 

= 1 + 2. = 1 + k\__-5— - 

Cq Cq x-^- + X,^ 

-° = ~ 1 _ k X X , 

C 2 p 



2 



which is equivalent to neglecting the terras with k, and 
k - tan cp sine i9- < 1 , or "better tan '-V< 1 . In addition; 



w + 0 r 

_o __2 _ _.o_ 

c' ""o ~ T 

I 0 ^ tan (p., - f,) ^ ; 



^0 - ;^_o ^^0 
2 
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or with 



^■0 - 



! 



4. ^ 0 
^ + 

° 2 



19. 
u^ 



— ° = (1 - -n.) X — (11) 



2 c 



0 riiX + \ 



iieve n, = induced efficiency of the hlade element at radi- 

'i ' 

us r of propeller in axial air flow of w^ velocity. 



Now equation (lO) hecovies; 
Ai =■ 



X + A ^ TliX + A lV X + a y 



i.p(i-^j) ,^.^>.^--\--i (10.) 

The determination of h is ootained frora Fi,^ure 4, 

c * / / c ^ 

where c_[_ is geometrically divided into — I- // — o. and 
2 2 2 

— S_ _L c__ "becatLse: 
2 2 



.1 



X _ _o ^ _i J? + 

2 2 2 

C ' C ' C ^ T 

When // and — i- ^ ~— we may put 

£i C> 2 O 



0 2 



c + 
0 2 



0 o 

2 2 



and 
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w ' 



cl Ai • _ o 2_ 

2 . c .: 



or 



w ' 



•' 0 



0 + 

2 



ci A i 4._ = 



Finally the last equations yield; 



Wo + 5l 



1.+ lL°- 

2_c_o 

Co 

1 + r- 



2 c, 



2c 



(12) 



Putting "b = 1 in equation (9) and Ai = ^ ^ 

I tn> equations (9) and (12) as first approximation from 
(10^), V7e can define..- b- "by means of the first two equa- 
tions. The thus oht'ained . h value inserted in (10 0 and 
the numerator and denominator of the latter developed ac- 
cording to k and ' Cl - ) , in which hoth quantities may 
"be cQxisidered as heing of the same order, then obtains to 



£\X zz — ^ • 

X -niX^H- p(i-'ni) Ax+A? 



X 



1 - kX 5^ - {l-r^O 

, X 4- A. n^:^ +x 



do") 

(the terms of the second and higher order "being neglected.) 

For the case of zero "load factor of blade element 
(or T^j^ = 1), its change in angle of flow becomes 



Aig = 



+ X" 



1 - k X 



X 1 



(10 -"a) 



(change in angle after neglecting t'He induced additive ve- 
locities, which may "be expressed as geometrical change in 
angl e . ) 
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But Ai drops considerably from this value \7ith the 
load rating (fig. 5) which indicates that the angular 
change is reduced "by the au t o r o t at i on . From this it is 
seen that highly loaded propellers are less sensitive to 
sideslip. Figure 5 shows, aside from the Ai curves de- 
fined from (lO'')» those accurately computed according to 
(8'), (9') and (12) at k = ± 0.2 for different \x as 
parameter versus T]^ , {\x ~ coefficient of advance at ra- 
dius x). It hecomes manifest that the exact Ai values 
are subordinate and "become greater only by higher load rat- 
ing. 

Now it should be of interest to examine more closely 
the distribution of the induced additive velocities due to 
sideslip u^^. To this end it is assumed that the angular 

change Ai is, in first approximation eqiial to the geomet- 
rical, i.e., to the neglected induced additive velocities 



Ai. = 



1 - k\ ^ ■ 

x^ + X^' 



Then, when expressing the corresponding velocity com- 
ponents as w*3 = w*- w^Q, u'g = u^ - u^Q 

and 



and approximately put: 



W ' = — ■* n ^ Tl ' =: X 



c' , u' - ^- 



and 

equation (9) yields 



\i + y.X E 



(9a) 
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The last equations give a clear picture of the distri- 
"bution of the induced additive velocities in the plane of 
the Dlade due to sideslip. Their ratio to the corresponding 
induced velocities of the propeller in axial air flow, i.e. 
when Uq^ = 0 and w = Wqi diminishes c ons ider ahly with 

the radius after having reached the maximum in the vicinity 
of X - \. The distribution over the periphery is in first 
approximation conformably to a sine curve. 

The total thrust and the torque are scarcely changed 
in their mean value hy the yawing position of the propeller, 
and can "be computed according to the usual method of the 
propeller theory. The two values increase in the half of 
the plane of the hlade where k = tan cp sine S- is positive 
and decrease in the other half, so that the mean value re- 
mains practically constant and equal to the propeller in 
axial air flow of Wq velocity and with the same revolu- 
tions. On the other hand, the center of gravity of the 
thrust is, under these circumstances, no longer in the pro- 
peller axis, as result of which a yawing moment perpendicu- 
lar to the axis of rotation is set up with respect to the 
propeller axis. Simultaneously there is a mean traiisverse 
force, perpendicular to the propeller axis, as a result of 
the unsymmet r ical distribution of the tangential forces over 
the periphery. 

The position of the resultant thrust which, moreover, 
is counted parallel to the propeller axis, is defined by 
coordinates of its intersection with the plane of the blade: 



/ / dS X R sin S- 




0 0 





X R cos d- 




0 0 




(13) 



The transverse force components are: 
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= /7 



dT cos S- 



2TT 1 

Ciz = / / dT sin 

0 0 



(14) 



In equations (l') and (2^, r\^^ ^ and T]^ ^ can be 
replaced "by their mean values and r]± and considered 

as i:eing equal to the respective quantities of the propeller 
in axial flow (constant over the periphery). It is, namely: 



a ^0 



a 



2 



wherein "b "becomes smaller or greater than 1, according 
to whether k is positive or negative and, for the rest, 
deviates only at large k values and in the neighborhood 
of the huh more appreciably from 1. This also holds for 

'Hi- 

Thus the replacement of Ai in (l') and (2') by 
equation (lOO reveals that 

f / dS X R cos i> = 0 

0 0 

and 

2TT 1 

//dTcosi5-=0 

0 0 

i.e., that the resultant of the thrust intersects the plane 
of the blade on the X axis (fig. l) and that the transverse 
force of the propeller lies on the Z axis. The vectors of 
the yawing moment and of the transverse force therefore coin- 
cide wi-th the Z axis. 

Developing the other two equations (13) and (l4), 
omitting the small terms of second and higher order and in- 
tegrating according to S- , then results in: 
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Yawing moment of propeller: 



° 2 0 X 0 r, 



4- ^ ^ 



X 2 + X= 



^ (2Tli - 1) - € (2Tli - 1) + T].\X - e\' 

+ p )?—^ 

(x ^ + /v^ ) (n^z 2 + p(i - ni)>.x + ) 



dx 



(15) 



wherein S - mean tlirnst, and the transverse force: 



X 
'Hi 



X ^ + X.^ ^ 



^ x^ + X^ 



+ -p X X, 



£ -1 (Srii - 1) x^ + (2ni - 1) x^ + e TiiXx + X 



2T 



(x^ +X^)(n,ix^ + p(l - T). ) X X 



dx 



(10 



To make the formulas clearer we insert 



x_ = ^ 



'Hi 



a 



^'x ^1 + X.^ 



c = M.__ 

R TTq R 



R Tiq 



and 



P — 



w 



2 

2 



cos^ cp 
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with ~ geometrical and X = ind-aced coefficient of 

advance on peripheral radius x and Crn = yawing moment 
and = transverse force coefficient. 

Therefrom follows: 



= 0 sin 2 cp 



(150 



wh ere 



(16') 



■5 



1 ^ ^'x 1 

4RTT i ^ \^ ^,^2 



1 + ^/T + X^3, y 1 + Xl^^ 



- £ + p 



2 ^K" ^'x 

(1 - e\'^) 



1 + p ( ^; 



xdx 



(15") 



2 + >; 



_ 1 



1 + £ — 



Hot e , that 



X' 



+ P 



and 



1 + p(\'x - X^) + x^'^ 
p = a - TT ' 



ax 



(IV) 



antu'l e 



of attack of profile in radians) are for infinite span. 

The last two equations offer a basis for the examina- 
tion of propeller energy and for its dynamic effect under 
practically any lihely running condition. 

With the present system of coor(3inates (figs. 1 and 2) 
a propeller rotating with velocity v/ at an angle to its 
axis and with angular velocity (i) is subjected to the follow- 
ing forces: a thrust in the Y axis (axis of rotation) con- 
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trary to axial component and a transverse force in the 

Z axis in direction of the transverse component u of 

^0 

the velocity of air flow, a torque about the Y axis con- 
trary to rotation go (and equal to the engine torque) and a 
yawing moment about the Z axis, "both - viewed contrary to 
velocity components w^ and u^^ - in the same effective 

direction. Thrust and torque are, as already shown, scarce- 
ly dependent upon cp when the axial component of Wq and oo 
remain constant, and equivalent for the projjeller in axial 
flow. These two principal forces are readily computed.* 

The transverse force like the yawing moment is depend- 
ent on the angle of attack of the propeller, in addition to 
being largely interdependent with all its parameters, that 
is, with the blade form and characteristics of the blade 
profiles, the coefficient of advance and loading and its 
course over the radius. These quantities being defined in 
the design for normal coefficient of advance, the desired 
forces are readily obtained by graphic integration from 
(15") and (16"). It is merely necessary to write the 

values in the respective formulas, which correspond to 
the propeller load curve corrected with K = f [r(1 - x) , 
z, sin)p*3» ^0 insure "finite number of blades." (Refer- 
ence 3.) For studies of propellers under operating atti- 
tudes at other coefficients of advance, (r,eference 4) 
Troller^s graphic method for propeller polars is of advan- 
ta.^e. He uses an auxiliary diagram (fig. 6) which shows 

c^ z t 

the -— — values for different \^ as parameter versus 
xR X 

angle j3 ' ^ - arc tan • basis of this diagram is 

the formula 

ia!.^ ^ 8tt 



^/l + (1 + XxX'x) 

from the general propeller theory, 
dc 

Rutting a = constant, which is acceptable for 

d i 

practically all profiles in the usual angle of attack range, 



*See reports cited here as well as others by Th . V. Bienen, 
Th. V. Xarman and Th. Troller. 
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zt 



and tracing corresponding to normal operating atti- 

tude in the graph, a straight line at slope 



d 








. xR J . 





d(3' di 



X 

c^ zt 

can "be drawn through this point, on which the — — — val- 

xR 

ue of all other operating attitudes lies. Even the exact 

[^Ca = ^(P*x^^ curves of the respective profiles can he 

traced for this purpose into the diagram in the — scale. 

xR 

Then we define through the intersection of the particular 

X line and the straight line ( or ^a^^.'^ curves the cor- 

V xR / 

responding and X'x value. To account for the "finite 

numher of "blades" the ordinates of the straight lines and 

c^z t 

— — curves are multiplied "by factor K - fC^^(l - ^) » z, 

sine |5 ' ' after which the and X^^^- values corres- 

ponding to the given Xx can "be read on the abscissa. In 
this manner the necessary X' , Co » / _ £_a ^ and e 

' Ca ; 

values are readily obtainable for any operating attitude. 

An example will illustrate the magnitude of the yawing 
moment and transverse force with a certain propeller as 
well as the shortest method of calculation: 

The dimensions of a propeller for a high-speed airplane 

ar e : 



Outside diam^eter 


D 




2 . 50 m 


Mean load factor 


a 




0.20 


Normal coefficient of 








advance 


X 




0 .312 


Produced thrust 


S 




2 90 kg 


At flight altitude 


H 




3,000 m 


With flight speed 


w 




2 90 km/h 


Engine power, sea level 


N 




5 50 hp 


Eevolut ions 


n 




1,980 r.p.m. 



15 



N.A.C.A. Technical Memorand-am Ho. 696 



The propeller was designed according to the Bienen- 
V. Karman method under the assumption of optimum thrust 
grading with a view to friction and finite "blade numoer. 
The "blade profile characteristics, constant across the 
whole radius, are: 

d c ^ 

Co = 0.80, „§: =4.9, p = — ^ =6.10, € = 0.02 

^ di ca 

TABLE I. 



X 


x= r/R 


.28 


.36 


.48 


.60 


.72 


■ .905 








Cazt(m) 


.394 


.400 


.400 


.359 


.288 


.158 










1.115 


.866 


.650 


. 520 


.433 


.345 










1.290 


.991 




. Dcsy 


/l D Q 


. OOO 






.312 


OX 














.063 


.0685 ^ 




.0295 


.0470 


.0775 


,100 


.111 


.090 








a y 
d X 


.0975 


.0950 


.0880 

1 1 


.0740 


.0560 


.0285 




\ 




c^zt(m) 


.180 


.175 


.180 


.170 


.153 


.0883 











1.43 


1.11 


.834 


.666 


.555 


.442 










1.53 


1.175 


.879 


.700 


.585 


.463 


! 


.40 


P 


13.10 13.95 


13.60 


12.95 


11.53 


10.95 


.043 


.066 




€ 


.037 


.040 


.039 


.037 


.033 


.031 








d <b 

u. 

dx 


.0196 


.031 


.0305 


.0672 


.0772 


.0620 






1 


d^i^ 
dx 


.0940 


.0910 


.0845 


.0720 


.0555 


.0280 


. 






Ca_zt(m) 


.564 


1 

.570 


.558 


.476 


.378 


1 

.210 








>^x 


.893 


.695 


.521 


.417 


.347 


.276 








>^'x 


1.107 


■ 

.854 


.630 


.500 


.421 


.328 






.25 


p 


3.34 


3.34 


3.43 


3.96 


4.0 


4.0 


.079 


.064 






.027 


.027 


.023 


.020 


.020 


.020 








dx 


.0365 


.0545 


.0945 


.123 


.139 


.117 








d^^' 
dx 


.0900 


.0885 


.0830 i .0705 


.0550 


.0290 
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Table I demonstrates a stiitaole calculation scheme. 
The calculation was also made for other operating attitudes 
( = 0.25 and \^ = 0.40). It is seen that whereas the 
coeTficient of the transverse force changes very little 
with the operating attitude, that of the yawing moment in- 
creases considerably at small coefficients of advance. 

For the normal X = 0,312 

= 0.063 sin 2 9 

and 

Cq = 0.0685 sin 2 cp 

or a yawing moment 

M = I? iiq c^ 39 kg/ m 
and a transverse force 

Q = E^TTq Cq = 34 kg 

for cp = 10^ and the prescribed operating,: conditions of the 
projjeller . 

The thrust is shifted along axis X "by the fraction 
of radius 

r 

= ~ = 13.5 per cent 

and the ratio of tre^nsverse force to thrust is 

11.7 -ner cent 
S ^ 

For comparison we have compiled in Tables II, III, and 17, 
and in Figures 7 and 8 the results of the experiments (ref- 
erence 5) made in the National Physical Laboratory, together 
with figu.res obtained from eq\iations (ibO and (16"). 
The model was a four-blade propeller with D = 0.305 m diam- 
eter and normal coefficient of advance X = 0.235 and a = 
0.365 load rating. Blade width, profile characteristics 
and coefficient of advance over the radius are given in 
Table II. 
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TABLE II. 



= r/R 



.15 



Gg^zt (m) 



P 
c 

d X 



.108 
1 . 568 
2.05 
3.50 

.060 

.077 



. 30 
.117 
.784 
.922 
4.75 
.042 

.097 



_^60_ 
.111 
.392 
.470 

6 . 50 
.032 

.086 



. 90 



.061 
.261 
.310 
6 . 95 
.025 

.0433 



.0725 



c^zt (m) 



X. 



P 
€ 

dx 



-4- 



.0690 
1 .00 
1 .102 
8 .05 
.056 
.093 



,0695 
. 50 
. 543 



.0411 
.333 
. 362 



10.55 I 10.30 

I 

.037 I .038 
.0437 



,071 



r 

.151 j 
.30 i 
. 4075 
4.77 
.030 

.0925 



1" 



CgZt (m) 



X 



X' 



e 

d X 



.164 
. 60 
.779 
3.39 
.059 

.1055 



.077 
.20 
.267 
5. 50 
.035 

.049 



.080 
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TABLE III. 



CP- 




5 


10 


15 


20 


2 5 


CalciJlat ed \ = 0 .235 




.0126 


.0248 


.0362 


.0465 


.0556 




X 


.244 


.22 6 


.222 


.217 


.213 


Measured 




.0136 


.0246 


.0356 


.0470 


.0587 



X 

Calculsted , 


TABLE 

.18 
.080 


.235 
.0725 


. 30 
.0710 




0 "X 

Measured at cp = 10 , 

Y 


r 

.150 |.171 
.0955.08 -0 


.187 
.0815 


.20 3 
.0755 


. 226 
.0720 


.247 
.0705 


.278 
.0370 


.292 
.0730 



The coinparison shows close agreeiLer.it iDetween our for- 
mulas and the experimental data. Tahle III and Figure 7 
show the Cq values for X = 0.235. Although the measure- 
ments did not exactly ohtain to this figure for \, it 
nevertheless was not necessary to correct the data, "because 
the coefficient of transverse force in this range is prac- 
tically unaffected hy \. (Fi^v- 8). When it is noted that 
our experiments were confined to snail k values (tancp< l) 
the amhit of agreement appears ahundant , the more so, since 
according to Figure 7, the discrepancies of the theoretical 
curve from the true attitude will not he abnormally large 
even ahove = 25^. 

The calculated and measiired values for variable 

coefficient of advance \ are compiled in Tahle lY. The ^|/' 
valu.es determined from (16") were arrived at from the data 
in Tahle II and Figure 6, whereas the test data were extra- 
polated conformahly to (16^). These data are plotted 
a^^ain-st \ in Figure 8 on the theoretically defined points. 

The load curve a - — - — evaluated from the measurements 
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is also given for comparison. 

It is manifest that the change in proportionality fac- 
tor i{/ of the transverse force for a given propeller of 
stated shape and dimensions with X is insignificant, as 
moreover, is readily seen from (15") and (16"). For 
certain propeller dimensions the lift coefficient of the 
hlade profiles tecomes lower "by increased X , and p = 
c ' 

a 

~~ correspondingly greater and vice versa. It follows 
from this, that when disregarding the inferior effect of 

the change in e and -Ijlx, the proportionality factor $ 

A, y 

of the yawing moment decreases as X increases, whereas 
that of ^\f remains practically constant. This fact is 
very significant when noting that the increase in X is 
followed "by a lower propeller-load factor (fig. 8) and hy 
a markedly higher ratio of transverse force to thrust. 

The very fact that c^ does not change with the co- 
efficient of advance and consequently with the flight 
speed, makes it possible to demonstrate the effect of the 
yawing propeller on the stability of the airplane very 
clearly. 

Assume, in place of the propeller, an airfoil of sym- 
metrical profile is rotatahly disposed about the propeller 
axis, that is, in such a way as to be able to automatically 
assume a setting perpendicular to the YZ plane (propeller 
axis - flight direction), and in the XY plane. Then the 
transverse force of the propeller can be replaced by the 
lift of this airfoil, provided this lift is always counted 
perpendicular to the propeller axis. For normal profiles 
in the usual range of angles, a fairly close approximation 
gives the lift coefficient as 

Ca - a sin 2 cp (a= constant) 

with cp = angle of attack, figured from the (cg^ - O) line. 
After approximating the usual geometrical profiles at a = 
2, the magnitude of this area (lift = transverse force) is 

f^ I t ^ TT = ^ Rtt 

^ a 

and the transverse force (lift): 
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Expressed in this form the term lends itself readily 
for use in the stability equations. 

Now in order to justify any general remarks about the 
behavior of a propeller in yaw the terms (15") and (16") 
must be rid of the propeller dimensions (zt and R) . 
This is accomplished by means of equation (l), which, 
transformed reads: 



or 



dS 



a - 



cazt 1 + \\ 1 - cX\ 



2rndr q 2 rrr 



X _ 



1 + \' 



2ttRx a — 



1 + 1 



eA. 



(la) 



The last equation, written in (15") and (l6"), 
leaves it dependent only on and and an 



(P = 



and C ) . 



Now it is readily apparent that the $ and values 
increase with \x for constant and profile character- 

istics (p = c'a and c). In other words, the ratio of the 
c^ and values to increases enormously with Xx 

under these circumstances. By constant \x and variable 
propeller loadings this ratio varies only slightly with the 
load rating and shows rather a tendency to drop at its in- 
crease. 



For the first approximation of the propeller force an 
old .er formula of the author can be used. To integrate 
equation (14) over the radius the blade number was assumed 
infinite and the total blade width expressed by 



zt = 8ttR ^ ^ 



as function of the radiua, which is ec^uivalent to the load 
dist ribiit ion 

a = 4 

^ 'He, ^i x^ 4- \^ 
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- -~ — = - re-placed eauation (11) and rio 

2 Co 2 Uo X - - 

and Tj^ were kept as constant across the radius. The ap-- 
proximation formula reads 



c = \lr a sin 2 .p 

where 



y = ^ In — — y 



2 



+ e n (2 ^ \ arc tan 4 +• 
'a V 



4- ___a ^ ^^^^ 1 

ca (1 + g-) 2 L 



I-±-4 ( ^+ S + CT].)]) (17) 

1 + X. 

and 

1 - 'Ha c»a 



s = 



The mean efficiencies can te determined from 

2 

V = J 

1 + V 1 + a 



r,i -Irrr—d - aX^ 

1 + v 1 + a 



with ^7 = mean load rating, 
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For the example in question the formula yields 

:= 0 .222 

or 

— = 0 .22 3 sin 2 9 
a 

For cp = 10^ angle of attack of propeller it is ^1 - 7.8 

a 

per cent against the more exact value of 11,7 per cent. 

RESUME 



The theoretical analysis reveals that a propeller in 
yaw is subject to axial thrust and torque and in addition 
to a transverse force and a yawing moment. The transverse 
force suhtends the axis in the plane defined "by it and the 
airflow direction, and acts in direction of the sideslip 
component "^q^ • "^^^^ force couj)le corresponding to the yawing 

moment lies in the plane through the propeller axis and a 
perpendicular to this and the direction of flow, and acts, 
when viewed opposite to the axial and transverse velocity 
component, in the inverse direction of the propeller rota- 
tion. For lightly loaded propellers with high coefficient 
of advance, that is, in high-speed airplanes, these forces 
can he comparatively great, and thereby affect appreciably 
the airplane stability and flight conditions during climb 
and banking. 

The transverse force may be assumed replaced by the 
lift of an airfoil disposed perpendicular to the flow di- 
rection of dimensions 1 t = JIL- s i n_2 j^-^ whose 

^a 

(cg^ ~ 0) line is coincident y/ith the propeller axis. 
This representation reveals very clearly the relation be- 
tween the transverse force of the propeller and the effect 
of the tail surfaces and their effect on the stability of 
the airplane. 



This effect is par t icrdar ly important in airplanes 
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with propellers mounted far ahead (or "back) of the center 
of £;ravity. The computed axial thrust and torque undergo 
no perceptihle change in mean value for the propeller in 
yaw, although its variation ahout the mean value which is 
of the order of magnitude of the transverse force and syn- 
chronized with the rotative speed, may induce dangerous 
"blade flutter as well as vibrations in the whole propeller- 
eng ine unit . 



II. EXFEEIIv'ENTS ON THE EFFECT OF THE SLIPSTREAM 
ON THE CONTROL SYSTEM OF AN AIRPLANE 



Part I of this report describes the phenomena produced 
in the plane of the "blade as result of the yaw of the pro- 
peller against the flow and relative flight direction. It 
was proved that there"by the air forces acting on the propel- 
ler can momentarily "be changed materially. 

Accordingly the flow pattern aft of the propeller will 
also "be different, so that the air forces on the parts of 
the airplane within the slipstream can he apprecia'bly af- 
fected. Ordinarily the mean direction of the developed 
slipstream, even apart from the spiral, deviates from the 
flow direction an-d alters, ahove all, the effective mean 
flow angle of the affected surfaces. But the flow picture 
is very complex and not very readily amena"ble to analysis. 

So in order to obtain a clear picture of the slipstream 

form and particularly of its influence on the controls, we 

explored the flow "behind the propeller and its effect on a 
fin hy experiment. 



1. Test Arrangement and Method 



The experiments were made on the propeller torque 
stand in the wind tunnel of the Aachen Aerodynamic Institute. 

The model propeller (D = 0.40 m diameter, fig. 9) was 
driven from an electric motor mounted on a post rotata"ble 
about two mutually perpendicular axes (figs. 10 and ll). 
The propeller axis was turned in the vertical plane against 
the flow direction. To minimize flow disturbances, we re- 
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sorted to fairin,--s. (See fis* 10.) The propeller revol\i- 
tioiis were recorded stroboscope. The thrust was recorded 
in axial direction "by a suitacle diaphragm or thrust meter- 
ing hoz. The instrument for recording the speed in magni- 
tude and direction is shown in ITigiires 12 and 13. The stat- 
ic pressure gauge A is rotatahle aoout two mutually per- 
pendicular, axes meeting in the test point, so that the po- 
sition of the test point remains unchanged. The opposite 
static tuhes of the instrument are connected to U tuhos 
B. When the axis of the static pressure ga^ige is coincident 
with the velocity direction the static pressures in one pair 
of tuhes equal that in the other pair. Then, when the U 
tuoes are set to zero, hy turning the instrument, the angles 
formed "by the air velocity with two mutually suht ended 
planes can "be read on dish C and plate D which are 
marhed off in degrees and vernier division. The zero posi- 
tion on "both degree divisions is determined in the same way 
by stopped propeller. The sensitivity of the test appara- 
tus for angle changes can "be regulated as desired hy using 
more or less sloping U tu'bes. For speed quantity measure- 
mxcnts the inner, shielded tube is connected with one of the 
outside to a manometer. 

For investigating the slipstream effect on a swing 
within its range, a fin with symmetrical profile of 600 mm 
span, 200 mm chord and 20 mm thickness was horizontally siis- 
pended on the test balance, so that the leading edge of the 
fin intersected the propeller axis at right angle, 1.10 m 
distant from the plane of the blade. The air loads on the 
fin were measured with suitable diaphragm boxes. 

The measurements were made for different angles of at- 
tach of the propeller axis against the air flow direction, 
i.e., 9 =^ 0^0, 5^, 10° and 15^, constant n ^ 2,500 r.p.m, 
and two air flow velocities corresponding to q = 7.5 and 
15 hg/m^ dynamic pressure, the axis being turned in the 
vertical plane. Magnitude and direction of the velocity in 
the slipstream was measured with the above described instru- 
ment in three successive cross sections parallel to the 
plane of rotation at z - 0.11, 0.58 and 1.16 m distance 
from the latter, that is, along one axis lying in the ver- 
tical plane and one perpendicular to this in the stream 
center. The instrument was carried in these two directions, 
as seen in Figure 14, and was adjustable parallel to the 
flow direction by means of a tube. The center of the 
stream could be accurately defined with the static pressure 
gauge, because of its quick reaction to the sudden direc- 
tional changes in velocity in this point. When the instru- 
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ment in the center itself is set parallel to the vortex axis, 
the manometer columns return to zero setting for reasons of 
symmetry and i^'ive very severe deflections even when only 
minutely displaced from its position. 

2. Results of Tests 



The primary ohject of the measurements was to gain an 
insight into the form of flow "behind a propeller. The re- 
sults are shown in Figures 15 and 16 for two running condi- 
tions. The velocities parallel to the stream center line 
are plotted against the distance from the stream center. 
The center line is defined "by the center points in the in- 
dividiial cross sections and shows at the same time the mean 
direction of the propeller slipstream. The slipstream re- 
mains, i]idependent from the angle of attack of the propeller, 
approximately cylindrical, and its "boundary with the dis- 
tance from the plane of the blade hecomes only little in- 
distinct "because of friction. The jet contraction is only 
sli>;'htly perceptible directly "behind the propeller, after 
which the slipstream slightly diverges again. 

The air speed in the plane of the "blade is the result- 
ant of inflow and induced velocity. As shown in Part I, 
the air forces of the "blade elements and as a result there- 
of the induced as well as the total velocities for the pro- 
peller in yaw are unsyrimet r i cally distributed over the pe- 
riphery. The phase of this distribution is shifted by the 
autorotation of the vortex system with the distance from 
the plane of the blade in relation to the radius and rate 
of slipstream rotation in the direction of rotatiQ^i of the 
propelier. This likewise alters the velocity distribution 
across the radius with respect to the distance from the 
plane of the blade. 

The mean direction of the slipstream is defined as fol- 
lows: the propeller experiences, aside from the mean axial 
thrust (leaving out the torque), a force component perpen- 
dicular to this, which in the present case is upward in the 
vertical plane. Accordingly the air quantity traversi};ig 
the plane of the blade is accelerated in both force direc- 
tions. The resultant mean air speed in the slipstream is 
therefore geometrically comiposed of the additive velocities 
corresponding to air flow velocity and the cited forces. 
Figure 17: the slipstream direction is graphically shown 
for both running conditions with the respective cp = 5^, 
10^ and 15^. 
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The transverse force for the model, computed according 
to the known method, was 

c = -2^- = 0.125 sin 2 cp 

^ R TTq 

The mean induced additive velocity corresponding to the 
thrust was taken from the measurements in cross section 
z - 0.11 m behind the propeller, wherehy 



for the other velocity component. 

For cp = 5° the thus defined mean slipstream direc- 
tion is in close agreement with the test data. With cp = 
10^ the measured angle between mean slipstream and air- 
flow direction is somewhat smaller, and with cp - 15^ it 
differs considerably from the calculated value. 

This discrepancy is due to the following: calculation 
did not allow for the fact that the propeller slipstream, 
when obliquely flowing into the main stream, is surrounded 
by the enveloping fluid, (Fig. 18.) Now there is a pres- 
sure difference on the area of discontinuity because of the 
probable vortex separation on the downstream side, which 
deflects the slipstream toward the air-flow direction. 
This effect is however, insignificant except at unduly 
large cp , so that in this range the mean slipstream direc- 
tion can be qiiite closely approximated by the above m.ethod. 

With the same workin^v conditions of the propeller we 
then investigated experimentally the effect of the slip- 
stream on the air forces directly on a fin with symmetrical 
profile (span; 600 mm, chord: 200 mm, thickness: 20 lam) , 
which, conformably to the horizontal tail unit, was placed 
horizontally about 1.10 m behind the plane of the blade, so 
that its leading edge met the propeller axis. 

Figure 19 shows the test data for the two operating 
attitudes with different 9 as parameter in the c^ - f(cCf) 
diagrams, Cg^ being the lift coefficient referred to total 
fin area and air-flow velocity, and af the angle of attack 
of the fin. When the propeller axis is parallel to the air- 
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flow direction the slipstream also is axial and the fin is 
in its center. Its angle of attack, hecaiise of the slixj- 
stream. rotation, increases on one side and. decreases on the 
other in the same measure. The result is an unsymmet r ical 
lift distr ihut ion (reference 6) ahout the longitudinal axis 
which is opposite to the propeller torque and "balances it 
in part. The mean lift however, remains ahout the same as 
if there were no rotation in the slipstream, and it may he 
considered as "being dependent on the axial velocity only in 
this case, i.e., that the fin loading; increases with the in- 
ducec axial velocity and the load rating of the propeller. 
But ordinarily, when the propeller axis slopes toward the 
flow direction, the mean rate of advance in the propeller 
slipstream also deviates from this direction, and thereby 
changes (that is, decreases) the effective mean angle of 
attack of the fin. On top of that, as Cp increases the 
fin emerges more and more from without the slipstream, so 
thaj^ in our case it is no longer in contact with it when 
cp =15^. (7igs. 15 and 16.) 

Based upon this reasoning the conditions may he brief- 
ly Slimmed up as follows: when the propeller axis is paral- 
lel to the direction of flow the lift of the fin is in- 
creased hy the increm.ent of the rate of advance in the pro- 
peller slipstream. Turning the propeller out of this posi- 
tion, the lift then deviates from this value, partly he- 
cause of the decreased effective mean angle of attack of 
the fin and partly as result of the emergence of the fin 
from within the range of the higher velocities. In Figure 
19 this is so displayed that the rise of the lift curves 
^a ~ f(af) decreases with Cp while simultaneously being 
shifted parallel to the axis of the abscissa. 

Another diagram traced by means of the (c^^ - f(a-p)] 

diagram in Pigure 19, is Figure 20, which shows the lift 
coefficient of the fin versus cp for different fin dis- 
placements with respect to the propeller axis af as 

0 

parameter. Here the rise of the [c^^ = f(?>)] is seen to 
increase with the load rating of the propeller, i,e., that 
the effect of the velocity increase in the propeller slip- 
stream: predominates. 

Any calculation of the air loads of a blade surface 
within the slipstream must first of all m.ake due allowance 
for r^nsymmet r ical velocity distribution and limited ranf::e 
of higher velocities, that is, of the same order of magni- 
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tude as the dimensions of the surface. In addition, the 
mean velocity far "behind the propeller is lower hecai^-se of 
the friction and its dist r ihut ion is altered. For parallel 
position of propeller axis to direction of flow, the ratio 
of the mean additive velocities at distance 0,11 m and. 1.16 
m "behind the plane of the "blade is - 1.25 for \ - 0.213 
and = 1.14 for \ = 0.300. The fact that the velocity dis- 
trihution far hehind is altogether different from that in 
the plane of the plate and load distr ihiat ion (see fi^'s. 15 
and 16) renders it, moreover, difficv.lt to find the inter- 
dependence between propeller loading; and control surfaces. 
Compiiting the lift of the fin according to 



A - c 



a / ^^^^P ^- 
0 ^ 



(13) 



where v/ - velocity far "behind the propeller, yields aonor- 
mally high figures. 

The ratio of mean velocity of advance in the slipstream 

far hejhind the plane of the Dlade to the undistiirhed flow 

velocity is by approximation according to the propeller 
theory: 



(a - mean load factor of propeller.) With the lift pro- 
port io.ial to the square of the speed and the measured value 

for a , equat ion 



yields 



— ( y . 

\w^ 



1 -f a 



(13) 



for A. - 



0.213 



0 . 30 - 



"a. 



2.05 
1 .42 



wherein c - lift coefficient of a supporting surface in 

"^0 

undisturbed flow and c_ in the sli-pstream, both referred 

cL - 



to the undisturbed flow velocity. 
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When extrapolating the c 
Figiire 19, which are referred 
600 X 200 mm, according to 



a 

to 



values from tlie plot in 
the total fin area I t 



(D = propeller diameter), the measurements yield with 



\ = I 0.213 " — 
^ 0 

0*30 - 



1.55 
1 .24 



Here the discrepancy "between the theoretical and the 
experimental c^^ values is markedly pronounced. It is 
largely due to the lower mean velocity in the slipstream as 
result of the friction, tecanse, when we substitute for w 
the mean velocity measured (on the fin) in the cross section 
1.16 m distant from the plane of the plate, and when we com- 



pute the corresponding 
have vvith 



values according to (18), 



we 



X 



0.213 



0.30 



1 . 53 
1 .20 



This shows that the air loads on a supporting surface 
within the slipstream are in close approximation proportion- 
al to the square of the mean velocity (referred to propel- 
ler disk area) at the particular point. However, this ve- 
locity is severely slowed up "by the air friction and very 
differently from the theoretical, defined according to the 
usual propeller theory, which makes mathematical analysis 
of the forces difficult. 

The effect of change of flow direction in the slip- 
stream is, as Figure 19 reveals, iviverse. 

Assum.ing the angle oetween flow and mean slipstream 
direction as proportional to 9 , that is: 



7 ^ ^ ^ 
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the 



effective mean angle of attack of the fin 



is : 



a 



=:cp-7-af = (1 



1^) cp 



a. 




o 



and "by constant velocity: 

Cq = constant, a = constant [(l = k) cp ~ ] 
a e 0 

The rise of the c^ (cp) curves is therefore diminished "by 
the slipstream deflection in the ratio of (l - k). The 
emer^:ence of the fin from without the slipstream acts iii the 
same sense. 

On the hasis of the last diagram (fic* ^O) it can 
finally "be ascertained that the staoility of the airplane 
is increased "by the effect of the free slipstream on the 
horizontal tail surfaces and that this effect grows with 
the load rating of the propeller. 



The quantity of air which the disk area described hy 
propeller "blades travels through, is accelerated rearward 
conforma"bly to the m.omentum of the forces acting on the pro- 
peller. This accelerated, air mass forms I'ehind the propel- 
ler the so-called slipstream., in which among others, the 
mean A^'elocity of advance is increased. 

If the propeller axis slopes toward the relative flight 
direction, the slipstream direction likewise changes, i.e., 
it sets up a downwash "behind the propeller. It was shown 
how this deviation can "be miathemat ically defined. The cal- 
culation is in close agreement with the experimental data 
except for abnormal slope of propeller axis to air flow di- 
rection; with appreciable sl07:)e the m.omentum of the slip- 
stream is partially split up in the s^ir rounding air "because 
of the circulation, and the downwash of the slipstreamx is 
d imin i shed . 

A control surface attached behind the propeller finds 
itself in the ran.^e of higher velocities and downwash by 
the yaw of the propeller, which increases with the slope. 
The effect of the velocity increase is much more pronounced 
than that of the downwash, as a result of which the dynamic 
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effect of a horizontal tail group in the slipstream is 
ther eby inc r eas ed . 

The deflection of the slipstream from the air-flow di- 
rection can, at least for small yav, he determined from the 
running conditions of the propeller. But the effect of 
speed increase cannot he resolved with the conventional 
methods of airfoil theory, because the flow is very differ- 
ent from the uniform potential flow. Besides, the velocity 
in the slipstream is not only slowed up hy the friction, 
hut its distrioution also is totally different at the place 
of the attached surface from that in the plane of the hlade 
and the thrust grading, as result of the coiling ujj of the 
spiral vortex surfaces. 



According to the formula 



A = c^^ / tdxp 

0 2 



Where w is, as before, the velocity far to the rear 
of the plane of the hlade, the lift of a surface placed in 
the slipstream is too high. To insure the exact amount due 
account must he taken of this quoted slowing up and differ- 
ent velocity dist r ihiit ion on the surface together with the 
fact that the range of the higher velocities is of the same 
order of magnitude as the dimensions of the control surface. 

The effect of the slipstream on the control surfaces 
and through it on the stability of the airplane can become 
comparatively great especially with highly loaded propellers, 
a fact which must be borne in mind in exact investigations 
of fli^^ht conditions and stability of airplanes. The exact 
determination of this effect is impossible until all the 
cited flow phenomena in the propeller slipstream are amena- 
ble to more precise theoretical and experimental analysis. 



Translation by J. Vanier, 
National Advisory Comm-ittee 
for Aeronautics. 
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Fig. 7 
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Fig. 15 Velocity dielrlbution in 
propeller slipstream. 
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Fig. 16 Velocity distribution in 
propeller slipstream. 
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